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Acoustic propagation boundary treatment error control. 
e Propagation by the Linearized Euler Equations (LEE). 
e Algorithms are 37@, 5% 7h and 9!” order in time and space. 


e Time accurate Damping Layer (DL) boundary treatment. 


Polynomial damping profiles are 2”¢, 4t2 6 and 8% order. 


Propagating solution uniformly set to zero on outer boundary. 


e Maximum relative absolute errors from O[10~%] to O[10~ 7]. 
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Linearized Euler Equation Conventions. 
e Cartesian coordinates (x,y) in two space dimensions. 
e Primitive variables V = (0,u,v,p)!, with o = 1/p. 
e Nondimensionalized with pr, Lr and speed of sound ap. 
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a 
© op =1/pr, tr =LpR/op, PR = 4ROR = a, and y= 1.4. 


e For both dimensional and nondimensional variables, 
a? = q = ypo, SO that the equation form is the same. 


Linearized Euler Equations (LEE) in 2D. 
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Vector form of the LEE in 2D with a source. 


AV; 
vA OV + Boe Vp aot : 


Vy. = —_— — (1/p,u,v,p)* , 


S(x,y,t) = (0,0,0,0.01 sin[2zt] exp[—25(x7n + y*)])/. 


e Either uniform base flow V, = (1,0.4,0,1/y)2, 
or a parallel jet V, = (1,0.4 + 0.4 exp[—25y2], 0, 1/4)". 


e The coefficient matrices A and B are from V~ 
while the coefficient matrices A, and B, are from Vj. 


Frame 001 | 24 WEWer4eh Wg | lee2dgsufc800dlczbd.f90 dat file 

LEE 2D: GS = 0.01 Sin[2 pi_ t] Exp[-25 (x+y 
c800 algorithm: h = 1/12, k = 1/36, ot = 7. 
DLC & ZBD for (x,y) in [-13,19] x [-18,18]. 


)], UF = (0.4,0). 


Pressure att = 25. 


p 
0.0005 
0.0001 
5SE-05 
1E-05 
5SE-06 

-SE-06 
-1E-05 
-SE-05 
-0.0001 
-0.0005 


Frame 001 | 24 May 2017 | lee2dgsgjo0diczbd.f90 dat file 
LEE 2D: GS = 0.01 Sin[2 Pi t] Exp[-25 (x+y’)], GJ = (0.4 + 0.4 Exp[-25 y’J,0). 


c800 algorithm: h = 1/12, k = 1/48 ot = 7. 
DLC & ZBD for (x,y) in [-13,24] x [-18,18]. Pressure att = 25. 


p 
0.0005 
0.0001 
5E-05 
1E-05 
5E-06 

-5E-06 
-1E-05 
-5E-05 
-0.0001 
-0.0005 


Damping Layers. 


e The numerical domain is Qn = [—x,,xpR] x [—yp, yr. 


Unsteady pressure source pg centered at the origin in &2yj,. 


Damping layers of width wp, wr and wy around Qj. 


The complete computational domain is 


Qc =[-(24,+ uz), ¢2R+uRMl x [-Wetvy)yrtwvyl), 
with Qa; C &24e as the inner core where accuracy is desired. 


Damping 


Dz(x) 


and in y, 


Dy(y) 


Nt? order polynomial damping profiles A. 


is done with polynomial damping profiles in z, 
= dp((a—-zpR)/we)®, for +ap<a<+antup, 
O, for —z;, <2 < +2p, 
= 6,((«@—az,)/wz)”, for —(a@,+wr)<«<—xz, 


= bp((y—yr)/wy)®, for tyrp<y<+yrtuy, 
QO, for — UB < Y < +YT, 
= 6p((y—yp)/wy)’, for —(ptwy) <y<—yp. 


Nt? order polynomial damping profiles B. 


We assume wR = wy = wy =w, Withhw=5 orw= 10. 


We assume dp = op = dy = 0, with O <6 < 50. 


We consider N=2, N=4, N=6and N=8. 


Damping is done with modified governing equations in Qo, 


© AV, = 


VA — +h V4 atts ee , + Pe (2) + Dy(y))V = 8S. 


Centered-Staggered /Cauchy-Kowaleskya/ Tailor 
cnoO Propagation Method. 


e Interpolation of an n2 data surface by dimensional recursion. 


e Cauchy-Kowaleskya recursion for time derivatives 


aerermay —— aeth+m(A,agV + BydyV) 
AxtOyeat™ AxOyat™ 
gethtm( Ad,V, + BayVp) 
7 Axtayrat™ 
OTT Dea Dy)V .. OF rors 
7 AxVOyat™ Axrayatm 


e Time advance by Taylor series. 
LO 


Maximum Relative Error with No Damping (6 = 0): A. 


e c400 computes with h=1/8 and k= 1/24. 
c6o0, c800 and c1000 compute with h=1/6 and k= 1/18. 


e Qn = [-3,7] x [-5,5], VY = (1, 0.4, 0, 1/7)". 


e Reported errors are 


_ max{|V(X) — V*(X)|: X in Quy} _ IV = V*llaw.ce 
max{|V(X)| : X in Qy} IV llan cc 


] 


ER ox 
where V* is a comparison solution from the same codes. 
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Maximum Relative Error with No Damping (6 = 0): B. 


Maximum relative error Erp. with w=5 and no damping. 


1.4930D-03 | 3.1436D-03 | 8.9365D-03 | 1.8767D-02 
50 | 6.6326D-05 | 2.7791D-04 | 7.3701D-03 | 2.4321D-02 


100 | 3.1044D-05 | 2.6507D-04 | 7.3527D-03 | 2.4521D-02 
200 | 2.3663D-05 | 2.2449D-04 | 7.2923D-03 | 2.4517D-02 
300 | 2.2813D-05 | 2.2534D-04 | 7.2911D-03 | 2.4516D-02 


Maximum relative error Ep. with w= 10 and no damping. 


50 | 5.2331D-04 | 9.1656D-04 | 2.3565D-03 | 1.0331D-02 
100 | 2.6468D-05 | 4.4113D-05 | 1.9533D-03 | 1.0577D-02 


200 | 1.9987D-06 | 1.2023D-05 | 1.9407D-03 | 1.0589D-02 
300 | 7.9275D-07 | 1.0818D-05 | 1.9408D-03 Te Re 
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Maximum Relative Error with No Damping (6 = 0): C. 


Outer boundary error is propagated like any other signal. 


ER oo increases with algorithm order, for each T’ and w, 
due to implicit damping from the greater diffusivity 
of the lower order algorithms. 


e ER ow generally decrease with simulation time T. 


A dissipitive dynamical systems with a periodic driving force 
converging to a periodic solution, any initial transient decays. 


LS 


Damping With a Uniform Flow: A. 


Uniform flow V, = (1, 0.4,0, 1/7)" for (x,y) in [—3, 7] x[—5, 5]. 


For c400, Ax = Ay=h=1/8 and At=k= 1/24. 


For c600, c800 and c1000, h= 1/6 and k= 1/18. 


Simulation time T = 100 or JT’ = 300 (or T’ = 200 for c1000). 


Damping layer widths w=5 or w= 10. 
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LEE2D GS UF(0.4,0) : c400, h = 1/8, k = 1/24, ot = 3, w= 5. 
Maximum Relative Error by Damping Profile Power N at T = 100 and T = 300. 
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LEE2D GS UF(0.4,0) : c600, h = 1/6, k = 1/18, ot = 5, w=5. 


Maximum Relative Error by Damping Profile Power N at T = 100 and T = 300. 


30 


Damping Amplitude delta 
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LEE2D GS UF(0.4,0) : c800, h = 1/6, k = 1/18, ot = 7, w=5. 


Maximum Relative Error by Damping Profile Power N at T = 100 and T = 300. 


Damping Amplitude delta 
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17 


LEE2D GS UF(0.4,0) : c1000, h = 1/6, k = 1/8, ot = 9, w=5. 


Maximum Relative Error by Damping Profile Power N at T = 100 and T = 200. 
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Damping Amplitude delta 
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Damping With a Uniform Flow: B. 


For N =2, Hao > 2) > Fpl: =0), 
and Ep(6 = 50) is O[102] greater than Ep(6 = 0). 


N =2 can produce errors larger than not damping at all. 

For N £2, Ep(6 > 2) is O[10] to O[102] less than Ep(é = 0). 
N polynomial damping has an Nt? order discontinuity at 0Q2j,. 
Quadratic damping more rapidly distorts the propagating so- 


lution surfaces as they enter the damping layer. 
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Damping With a Uniform Flow: C. 


Smallest observed error Ep with w=5 at T’= 100. 


algorithm | _N=2 | N=4 | N=6 | N=8 


2.1854D-05 | 6.4663D-06 | 6.5709D-06 | 6.2108D-06 


3.1538D-05 | 6.5983D-06 | 6.4781D-06 | 7.0472D-06 


5.4334D-05 | 8.2486D-06 | 6.3638D-06 | 7.2109D-06 
1.3842D-04 | 1.0321D-05 | 5.9463D-06 | 7.4360D-06 


Smallest observed error Ep with w=5 at T= 300. 


algorithm | N=2 | N=4 | N=6 | N=8 | 
1.6771D-05 | 8.9120D-07 | 7.2492D-07 | 7.2492D-07 
2.6601D-05 | 9.9872D-07 | 7.3028D-07 | 8.1127D-07 


4.8991D-05 | 2.7558D-06 | 7.4484D-07 | 1.3917D-06 
1.3421D-04 | 3.2921D-06 | 1.7333D-06 | 2.6527D-06 
C1000" at. = 200. 
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Damping With a Uniform Flow: D. 


O[10-*] < Ep(T = 100) < O[10~°], 
O[10—-4] < Ep(T = 300) < O[10~J. 


The upper right section has the lowest levels in each table. 


Higher N works better with all of the algorithms, 
with more necessary improvement for higher order algorithms. 


As a rule of thumb, the damping profile power should be at 
least as large as the algorithm order. 


The 9£ order c1000 might need either N > 10 or T > 300. 
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The Effect of Damping Layer Width: A. 


Uniform flow V, = (1, 0.4,0, 1/7)” for (z, y) in [—3, 7] x[—5, 5]. 


c800 algorithm with h= 1/8 and k = 1/24. 


Damping layer widths w=5 and w= 10. 


Simulation times T’= 100 and JT’ = 300. 
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LEE2D GS UF(0.4,0) : c800, h = 1/6, k = 1/8, ot = 7, w=5 and w= 10. 


Maximum Relative Error by Damping Profile Power N at T = 100. 


Damping Amplitude delta 


——— 


wo-Oee 


N2T100 w05 
N4T100 w05 
N6T100 w05 
N8T100 w05 
N2T100 w10 
N4T100 w10 
N6T100 w10 
N8T100 w10 


20 


LEE2D GS UF(0.4,0) : c800, h = 1/6, k = 1/18, ot = 7, w = 5 and w = 10. 


Maximum Relative Error by Damping Profile Power N at T = 300. 
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The Effect of Damping Layer Width: B. 


Error reduction by doubling w is greater for lower JN. 


Error reduction by tripling Z’ is greater for higher JN. 


For N = 2, error decrease is more efficient by increasing w, 
but the error is generally the largest if all else is the same. 


For N #2, error decrease is more efficient by increasing T’, 
and most efficient by increasing both w and T. 
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Damping With a Parallel Jet: A. 
V, = (1, 0.4+).4 exp[—25y7], 0, 1/7)? OR Vy, = (1, 0.4, 0, 1/4). 
Qn = [-3,9] x [—-5,5] OR Qy = [-3, 7] x [—5, 5]. 
For c800 with h= 1/12 and k= 1/36. 
Simulation time 7’ = 50. 


Damping layer width w = 5. 
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LEE2D GS UF(0.4,0) and GJ(0.4+0.4 Exp[-25y“2],0), c800, h = 1/12, k = 1/36, ot = 7. 


Maximum Relative Error by Damping Profile Power N at T = 50 with w = 5. 


Damping Amplitude delta 


N2T050uf 
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Damping With a Parallel Jet: B. 


This data is for the relatively short simulation time T’ = 50. 


This grid resolution is doubled from the uniform flow simu- 
lations. 


Errors from the uniform and jet flow differ by + 2 to #5, 
so the boundary treatment performs similarly for both. 


All of the effects of damping profile power, damping layer 
width, simulation time and algorithm order should apply to 
variable coefficient jet simulations with these algorithms just 
as they did to constant coefficient uniform flow simulations. 
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Conclusions. 


Error can decreases as simulation time T' increases, 
from decaying transients. 


Error tends to decrease as the damping layer width increases, 
more effectively and efficiently for smaller N. 


Error tends to decrease as damping power JN increases, 
more effectively and efficiently for larger N. 


As a rule of thumb, the damping profile power N 
Should be at least as large as the algorithm order. 
29 


Quadratic damping produces the largest errors. 


Damping with a jet and a uniform base flow are similar. 


Dissipitive algorithms implicitly dampen. 


High order algorithms propagate accurately, 
both intended signals and also errors. 
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